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In this paper two tlieoretical approaclies for tlie calculation of the rate of quasi-stationary, two- 
dimensional magnetic reconnection with nonuniform anomalous resistivity are considered in the 
framework of incompressible magnetohydrodynamics (MHD). In the first, "global" equations ap- 
proach the MHD equations are approximately solved for a whole reconnection layer, including the 
upstream and downstream regions and the layer center. In the second, "local" equations approach 
the equations are solved across the reconnection layer, including only the upstream region and the 
layer center. Both approaches give the same approximate answer for the reconnection rate. Our 
theoretical luodel is in agreement with the results of recent simulations of reconnection with spatially 
nonuniform resistivity by Baty, Priest and Forbes (2006), contrary to their conclusions. 
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I. INTRODUCTION 

Magnetic reconnection is one of the most important 
processes of plasma physics, and is believed to play the 
central role in observed phenomena in laboratory and 
cosmic plasmas. At the same time there has been a long- 
standing debate about the correct theoretical model of 
magnetic reconnection. Most previous theoretical and 
numerical work focused on reconnection processes in two- 
dimensions, in which all physical scalars and vectors are 
independent of the third coordinate (z). There exist two 
original and well known models of magnetic reconnec- 
tion with constant resistivity. First, the Sweet-Parker 
reconnection model [l|, which predicts a slow mag- 
netic reconnection rate in hot low-density plasmas. Sec- 
ond, the Petschek model in which a fast reconnec- 
tion rate is achieved by introducing switch-off magneto- 
hydrodynamic (MHD) shocks attached to the ends of a 
relatively short reconnection layer in the downstream re- 
gions. Many numerical simulations have been carried out 
to discriminate between these two models. More recent 
high-resolution simulations generally favor the Sweet- 
Parker model of slow reconnection in the case of con- 
stant resistivity, and do not confirm the Petschek theo- 
retical picture for the geometry of the reconnection layer 
with shocks [1, H, @. However, reconnection becomes 
much faster and Petschek-like if resistivity is not con- 
stant and is enhanced locally in the reconnection layer 
(numerical studies of this case were pioneered by Ugai 
and Tsuda 0, [1], by Hayashi and Sato d, [13], and by 
Scholer 

These results have been called into question in a recent 
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paper, Baty, Priest and Forbes (2006) (l5|, which claims 
that Petschek reconnection can occur when the resistiv- 
ity is not absolutely constant but varies by an arbitrarily 
small amount. This paper has motivated us show on the 
basis of an earlier paper of ours, Malyshkin et al. p^ . 
that the Petschek shocks are produced not by the varia- 
tion of the resistivity but by the rate of variation which 
for their case is very large and unlikely to happen natu- 
rally. 

In 2001 Kulsrud [l^l provided some theoretical insight 
into magnetic reconnection process which qualitatively 
explained the results of recent simulations. He considered 
quasi-stationary, two-dimensional magnetic reconnection 
in the classical Sweet-Parker- Petschek reconnection layer 
with zero guide field (B^ = 0) , zero plasma viscosity and 
an anomalous resistivity that is a piecewise linear func- 
tion of the electric current. Kulsrud was first to suggest 
that one has to calculate the half-length of the recon- 
nection layer L' from the MHD equations and the jump 
conditions on the Petschek shocks, instead of treating L' 
as a free parameter (as Petschek erroneously did when 
chose L' to be equal to its minimal possible value under 
the condition of no significant disruption to the plasma 
flow). As a result, in the case of constant resistivity Kul- 
srud correctly estimated the layer half-length L' to be 
approximately equal to the global magnetic field scale, 
L' K, L. In this case the Petschek reconnection rate re- 
duces to the slow Sweet-Parker reconnection rate [l^, [13 , 
in agreement with numerical simulations. The second re- 
sult obtained by Kulsrud is that in the case when resistiv- 
ity is non-constant but anomalous and enhanced (e.g. by 
plasma instabilities), the reconnection rate can become 
considerably faster than the Sweet-Parker rate. 

In our recent paper (Malyshkin, Linde and Kul- 
srud [13])) to which we will hereafter refer as MLK2G05 
paper, we put Kulsrud's derivations on a rigorous analyt- 
ical basis and extended his model by using a new theo- 
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retical approach to calculate the reconnection rate. This 
approach is based on "local" analytical derivations across 
a thin reconnection layer, and it is applicable to the case 
when resistivity is anomalous, nonuniform and an arbi- 
trary function of the electric current and of the spatial 
coordinates. We included the case of non-zero guide field 
Bz and non-zero plasma viscosity in our model. We 
found an approximate formula for the reconnection rate 
which confirmed Kulsrud's theoretical results. 

The present paper has two goals. First, to calculate 
the reconnection rate, and second, to compare it with 
the recent simulations of Baty et al. [15]. 

In the next section we simultaneously follow two theo- 
retical approaches to the calculation of the reconnection 
rate. In the first approach, which we call the "global" 
equations approach, we derive and solve approximate 
MHD equations for a whole reconnection layer, includ- 
ing the upstream and downstream regions and the layer 
center. These theoretical derivations are similar to those 
done before, except for an important difference. Namely, 
we take into consideration an additional important equa- 
tion, the spatial homogeneity of the z-component of elec- 
tric field Ez along the reconnection layer. This equation, 
together with the jump condition on Petschek shocks, al- 
lows us to find the reconnection layer length L', which 
must be determined consistently |12l |. The second the- 
oretical approach, which we call the "local" equations 
approach, basically coincides with the calculations done 
in the MLK2005 paper. In this approach we derive and 
solve approximate MHD equations across a thin recon- 
nection layer, including only the upstream region and 
the layer center. Both approaches give the same approx- 
imate formula for the reconnection rate, which is valid in 
the general case of an arbitrary nonuniform anomalous 
resistivity [see Eq. pT|) ]. Using these two approaches si- 
multaneously allows us to better understand the physics 
of a reconnection process. 

In Sec. |TTT] we pursue the second goal of this paper. 
We consider the case treated by Baty et al. [l^, when 
resistivity is a prescribed function of the two spatial 
coordinates x and y (remember that we consider two- 
dimensional reconnection, so that no physical quantities 
depend on z). We argue that our theoretical model is 
in agreement with the results of these recent simulations 
of reconnection by Baty, Priest and Forbes [iB]. This 
agreement contradicts their theoretical conclusion that 
Petschek shocks exist with constant resistivity, and we 
explain why. Finally, in Sec. IIVI we discuss our results. 



II. RECONNECTION WITH NONUNIFORM 
ANOMALOUS RESISTIVITY 

In this section we consider magnetic reconnection with 
nonuniform anomalous resistivity in the classical two- 
dimensional Sweet-Parker-Petschek reconnection layer, 
shown in Fig. [TJ The layer is in the x-y plane with the 
X- and y-axes being perpendicular to and along the layer 
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FIG. 1: (Color online) The geometrical configuration of the 
classical Sweet-Parker-Petschek reconnection layer is shown. 
Petschek shocks exist only if the reconnection is considerably 
faster than the Sweet-Parker reconnection rate. 



respectively. The length of the layer is 2L', which is ap- 
proximately equal to or smaller than the global magnetic 
field scale L that will be introduced below. The thick- 
ness of the layer, 25 o, is much smaller than its length, 
i.e. 25o ^ 2L'. The classical Sweet-Parker-Petschek re- 
connection layer is assumed to have a point symmetry 
with respect to its geometric center point O in Fig. [1] 
and reflection symmetries with respect to the axes x and 
y. Thus, for example, the x- and y-components of the 
plasma velocities V and of the magnetic field B have the 
following simple symmetries: 14(±a;,=Fy) — ^yxix^y), 
Vy{±x,Ty) = TVy{x,y), B^{±x,Ty) = TBx{x,y) and 
By{zLx,^y) — ±By{x,y). There could be a pair of 
Petschek shocks attached to each of the two reconnec- 
tion layer ends in the downstream regions. Because of 
the MHD jump conditions on the Petschek shocks, there 
must be a nonzero perpendicular magnetic field B' in 
the downstream region at point O' in Fig. [1] [H, [13. 
If the plasma viscosity is small, then the plasma out- 
flow velocity Vq^j in the downstream region at point O' 
is approximately equal to the Alfven velocity Va calcu- 
lated in the upstream region at point M (refer to Fig.[T]). 
The plasma inflow velocity Vr in the upstream region at 
point M , outside the reconnection layer, is much smaller 
than the outflow velocity, Vr <C Va- Finally, the mag- 
netic field Bm at point M, outside the layer, is mostly in 
the direction of the layer (i.e. in the y-direction). 

Now let us list four assumptions that we make about 
reconnection process. First, we assume that the charac- 
teristic Lundquist number is large, which (by our defini- 
tion) means that resistivity is negligible outside the re- 
connection layer. Second, we assume the plasma flow is 
incompressible. Third, for simplicity we neglect plasma 
viscosity. (The case of nonzero viscosity is treated in 
the MLK2005 paper Fourth, we assume the re- 

connection process is quasi-stationary. This is true if 
the reconnection is slow, Vr/Va ^ 1, and that there 
are no plasma instabilities in the reconnection layer. For 
an incompressible viscousless plasma the assumption of 
slow reconnection is equivalent to the assumption that 
the reconnection layer is thin, 5o/L' « Vr/Va ^ 1. The 
above assumptions are standard in the Sweet-Parker and 
Petschek reconnection models. Please note that we make 
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no assumptions about the values of the guide field Bz . In 
other words, our derivations apply to what is called "two- 
and-a-half dimensional" reconnection. 

For brevity, we use physical units in which the speed 
of light and four times tt are replaced by unity, c — 1 and 
Att = 1 In these units the MHD equations that we 
need to find the reconnection rate are as follows. Fara- 
day's law V X E = —d'B/dt for the x- and y-components 
of magnetic field in two dimensions is 

dEz/dy = -dB^/dtKiO, dE^/dx = dBy/dt ^ 0, (1) 

where 9B/9t « because of the quasi-stationarity of 
reconnection. From Eqs. ([T]) we see that the electric field 
z-component is constant in space, i.e. = Ez{t) is a 
function of time only. Next, neglecting the displacement 
current in the framework of MHD, Ampere's law for the 
z-component of the current is 

jz = (V X B), = dBy/dx - dB^/dy. (2) 

Ohm's law for the spatially uniform z-component of the 



electric field is 

Ez{t) — —VxBy + VyBx + rjjz — constant in space, (3) 

where resistivity rj — r]{jz,x,y) is an arbitrary func- 
tion of the electric current z-component and of the two- 
dimensional coordinates p^ . Next, the equation for 
plasma acceleration along the layer in the y-direction is 

p(y ■ V)Vy = ~{d/dy)[P + Bl/2 + Bl/2] + (BV)B„(4) 

where p is the constant plasma density, P is the sum of 
the plasma pressure and the guide field pressure Bl/2. 
In addition, we have 

dxVx + dyVy = 0, dxBx + dyBy = (S) 

because the field and the velocity are divergence- free. 

Now, from MHD equations (H))-® we obtain the fol- 
lowing table of "global" and "local" equations for the 
reconnection layer: 





Global Equations 


Local Equations 


1**', Ampere's law 


jo = {jz)o « idxBy)o ^ Bm/So 


2"^^, Incompressibility 


VrL' « V^^.So 


{dyVy)o = -{dxVx)o « Vr/6o 


3'^'^, Plasma acceleration 




p{dyVy)o^ = -{dlP)o+3o{dyBx)o 


4**^, Jump condition on shocks 






5"^, Ez = const across the layer 


Vojo = VRBrn 


6*'', Ez = const along the layer 


Vojo = i]'z + K'ut^i 


Q = dl{m.), + 2{dyVyUdyBx)o 


7*'', Unknown quantities 


jo, 6o, Vr, i', FJut, B'^ 


jo, 5o, Vr, {dyVy)o, {dyBx)o 



(6) 



Here quantities with subscript o are taken at the re- 
connection layer central point O, while quantities with 
the prime sign are taken at point O' in the downflow 
region (refer to Fig. [T]). The column "Global Equa- 
tions" includes equations that are written at points O, M 
and O', these equations represent the "global" equations 
approach to the calculation of the reconnection rate. The 
column "Local Equations" includes equations that are 
written only at points O and M and represent the "lo- 
cal" equations approach. Note that equations on lines 1 
and 5 enter both columns in the same form. The first 
line of the table includes the Ampere's law equation ([2]), 
with the dyBx term neglected because it is small, and the 
dxBy term estimated at point O. The second line of the 
table includes the plasma incompressibility condition ([5|) 
in its "global" and "local" forms. The "global" form 
is the mass conservation equation, while in the "local" 
form the dxVx term is estimated at point O. The third 
line contains equations for plasma acceleration given by 
Eq. (gl). The "global" equation V^^^ « Va reflects the 
well known result that in the absence of viscosity the 



plasma outflow velocity in the downstream region is ap- 
proximately equal to the Alfven velocity calculated in 
the upstream region [J, 0, 0| • The "local" equation for 
plasma acceleration results from differentiation of Eq. (j4|) 
with respect to y and taking the symmetries of the recon- 
nection layer into account. The jo{dyBx)o term in this 
equation is the magnetic tension force, while the pressure 
term is equal to 

{dlP)o = dl{Bl/2)^ = B^idlByU < 0. (7) 

Thus, the drop of pressure P (which includes the plasma 
and guide field pressure) along the layer is equal to the 
pressure drop of the parallel magnetic field outside the 
layer. This result follows from the force balance condi- 
tion for the plasma across the reconnection layer (in anal- 
ogy with the Sweet-Parker derivations), and its rigorous 
proof can be found in the appendix A of the MLK2005 
paper [ij]. The last equality in Eq. ([7]) comes from the 
layer reflection symmetry with respect to the x axis. 

Next, the fourth line in Table [5] includes the standard 
jump condition on the switch-off MHD Petschek shocks 
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attached to the ends of the reconnection layer [l^, ■ 
This condition is a "global" equation, and the plasma 
incoming velocity at point M' is estimated as being 
approximately equal to the plasma incoming velocity Vr 
at point M (refer to Fig. [T]). There is no corresponding 
"local" equation because we do not consider the down- 
stream region and Petschek shocks in the "local" equa- 
tions approach! Equations in the fifth and sixth lines 
of Table [6] directly result from Ohm's law equation ([3]). 
Namely, we use the spatial homogeneity of Ez (t) , which is 
a consequence of the quasi-stationarity of reconnection. 
To obtain the single equation in the fifth line and the 
"global" equation in the sixth line, we equate the Ohm's 
law expression for at points O, M and O' . To obtain 
the "local" equation in the sixth line, we take the second 
order partial derivative dy of Eq. ([3]) at point O. Finally, 
the seventh line in Table [6] lists all unknown physical 
quantities to be estimated in the "global" and "local" 
equations approaches. Note that quantities jo, So, Vr, 
{dyVy)o and [dyBx)o are "local" (i.e. defined at the layer 
central point O and at point M in the upstream region), 
while quantities L' , V^^^ and B'^ are "global" (i.e. defined 
at point O' in the downstream region) . 

There are a few additional equations that we need. 
First, we use the second-order Taylor expansion of -qjz 
along the y-axis to estimate the term in the 6th line 
in Table [g 



connection layer central point O: 



VoJo + (i'V2) {jo{dlv)o + 
+ [Vo + 3o{dr]/djz)o]{dyjz)o} 



(8) 



Note that the first-order Taylor expansion terms are zero 
because of the symmetry. Second, we define the field 
global scale L and the resistivity scale as (20j 



/{d^yByU, ll = ^2TJo/{d'yV)o. (9) 



Third, the y-scale of the current jz, to a factor of order 
unity, turns out to be about the same as the y-scale of 
the outside magnetic field. 



3o\dljz)o^B~\d^yByU 



(10) 



This result can be understood by taking the second order 
partial derivative dy of the Ampere's law equation jo ~ 
Bm/So given in the first line in Table [6l while keeping 
5o constant because the partial derivative in y is to be 
taken at a constant value oi x — So- The detailed proof 
of Eq. Uni) is given in the appendix B of the MLK2005 
paper. 

Now we have all equations necessary to find the re- 
connection rate and all other unknown physical param- 
eters. In the "Global Equations" column in Table [6] 
we have six equations and six unknowns, and in the 
"Local Equations" column we have five equations and 
five unknowns. Using the "local" equations and equa- 
tions jll), dH) and (Uni), we obtain the following approx- 
imate algebraic equation for the z-current jo at the re- 



3 + :^ 



drj 
djz 



hi 



A m 



(11) 



where resistivity rjo = ri{jz — jo, x — 0,y — 0) is a func- 
tion of jo. The "global" equations give a similar result 
with 3 replaced by 1 in Eq. pT|) , which is a less accurate 
result due to additional approximations made in Eq. ([8]) 
and in the formula « Vr in the fourth line of Table [S] 
Given the resistivity function = rj{jz,x,y), as well 
as the magnetic field Bm and its global scale L = 
—2Bra/ {dyBy)m both calculated at point M in the up- 
stream region, we can solve the algebraic Eq. pT|) for 
the current jo- Once jo is calculated, we can easily find 
the reconnection rate, which is the rate of destruction of 
magnetic flux at point O and is equal to the electric field 
z-component Ez{t) = rjojo- We can also find all the other 
physical parameters, by using the equations in Table [SI 



So 

Vr 

idyVy)o 
{dyBj;)o 



Bm I jo, 

Vojo/Bjn < Va, 
Vofo/Bl « Vr/6o, 

Jo iylJVl ~ 2BI/jIL^) 



L' « So 



Va, 
Va 



BmiVR/VA), 



Va 



BL 



Vr {dyVy)o {dyB.,)o 



(12) 

(13) 
(14) 

(15) 
(16) 
(17) 

(18) 



The last approximate equality in Eq. psp is valid be- 
cause the second term on the right-hand-side of Eq. (fTS]) 
can be neglected. This is because 'iB'^J j^L"^ < V^/V^, 
which follows from Eqs. (fT5|) and (fTT|) (also refer to foot- 
note [28]). 

Equations (fTT|) - (fT8l) are very general results for quasi- 
stationary magnetic reconnection with no assumptions 
about the functional form of resistivity and the guide 
field value. Now let us look at the three terms on the 
left-hand side of Eq. pT|) . When resistivity is constant, 
the only term left is the first term. The second term is 
clearly related to the dependence of anomalous resistivity 
on the current. The third term becomes important when 
resistivity is ad hoc localized in space. As a result, in 
the end of this section we consider three special cases of 
magnetic reconnection, in which Eqs. pT|) - p^ reduce to 
simpler formulas. These three cases correspond to dom- 
ination of the first, second and third terms respectively 
on the left-hand-side of Eq. ifTTj) . and they are as follows. 

The first case is the quasi-uniform resistivity case, 

if !■> ma,x[{jo/r]o){d'n/djz)o,L^/ll], then 

jo~{Bm/L)S^'^ where S=VALho, (19) 
Vr/Va^S-^'\ 5o^LS-^I\ L'^L. 

Here we introduce the Lundquist number So = VaL/tIo 
and assume for our estimates that 3^/^ ~ 1. Equa- 
tions are the familiar Sweet-Parker results |ll, 0]- 



5 



Thus, if resistivity is quasi-uniform or uniform, then re- 
connection is Sweet-Parker. 

The second case is the Petschek-Kulsrud reconnection, 

if Oo/»7o)(9r//9j^)o > niax[l,L2/;2]^ then 

VrIVa ~ ^olL' « {{B,jYAL^){d^ld3.)o\ , (20) 

This is the case of fast reconnection with Petschek ge- 
ometry and shocks. Equations ([20|) were first derived 
by Kulsrud [l^l- This is the reason we caU this case 
Petschek-Kulsrud reconnection. 

The third case is the case of reconnection with spatially 
localized resistivity, 

if L'^/ll:$> mSix[l,{jo/r]o){dri/djz)o], then 

jo ~ {B„Jlr,)Sl^^, where S^j = VaIti/iIo, (21) 

Vr/Va « s;^^'\ So « ir,s;,^'\ l' « « l. 

Here we introduce the effective Lundquist number 5*^ = 
VaIti/Vo that is based on the resistivity scale Irj given by 
Eq. ^ . Equations (PT|) are the same as the Sweet-Parker 
equations with the global field scale L replaced by 
the resistivity scale Irj. When resistivity is strongly lo- 
calized, Irj <^ L, the reconnection becomes fast yielding 
Petschek geometry and shocks. 

We postpone the detailed analysis and discussion of 
our theoretical results until the last section of the paper. 

III. RECONNECTION WITH SPATIALLY 
LOCALIZED RESISTIVITY 

In this section we compare our theoretical results to 
the results of recent simulations of reconnection with spa- 
tially nonuniform resistivity by Baty, Priest and Forbes 
(2006) They took the resistivity as 

Tjix, y)^m + iVo - m) exp [-{x/Q^ - {yllyf\ (22) 

and considered different values of parameters r^o, ryi, 
and ly. They found that the reconnection rate does not 
depend on the value of Ix- This is in agreement with 
our Eq. (|lip . which includes only dyr]{x, y) derivative via 
resistivity scale Ij^ given by Eq. ®- As for the depen- 
dence on the other parameters, in all their simulation 
runs Baty et al. found the Petschek geometrical configu- 
ration with shocks and reconnection rate faster than the 
Sweet-Parker rate. They paid special attention to the 
case when r]o — Vi Eq. (|22p is small and the resistivity 
is weakly nonuniform. They called this case a "quasi- 
uniform resistivity" case, and observed the Petschek so- 
lution in this case as well. At the same time, we theo- 
retically derived the Sweet-Parker solution for the quasi- 
uniform resistivity case, given by Eqs. p^ . Thus, there 
is a disagreement between our theoretical results and the 
claims by Baty et at, which needs to be addressed. 

The reason for this disagreement is that the resistivity 
used by Baty et al. was not actually quasi-uniform in all 



their simulation runs. In fact, when resistivity depends 
only on coordinates, the condition of a quasi-uniform re- 
sistivity is 1 3> L'^/l'^, refer to Eqs. (fT9|) . In other words, 
the resistivity localization scale Ir/ must be much larger 
than the field global scale L. In their most uniform 
resistivity simulation run Baty et al. used rjo = 10"'*, 
771 = 9.3x10"^, ly = 0.1 and L « 1. According to Eq. 
in this case Irj = ly\/Vo/ iVo ^ Vi) = 0.378, which is con- 
siderably smaller than L w 1. Thus, while the resistivity 
function rjix, y) = 9.3 x lO^^^ _,_ 7 ^ iQ-e exp[-(y/0.1)2 - 
{x/l^f'] would indeed appear close to uniform when its 
graph is plotted, its second derivative dyii is large due 
to small value of ly = 0.1. As a result, this resistiv- 
ity should be viewed as far from uniform and as rather 
well localized, /l'^ « 7. In their other simulation runs 
Baty et al. used even stronger localized resistivity with 
larger values of rjo — 771 in Eq. ((22|) . In fact, they were not 
able to run any simulations with strictly uniform resistiv- 
ity 77 — const because of an instability of the reconnection 
layer. The reason for this instability possibly lies in the 
boundary conditions that Baty et al. used, which might 
conflict with our Eq. (fTT|) . The latter must be satisfied 
for quasi-stationary reconnection. 

We would like to quantitatively compare our theo- 
retical results to the results obtained by Baty et al. in 
their most uniform resistivity simulation run. They had 
l„ = 0.378, L w 1, 770 = 10-4, Va « 1, « 0.9 
and they found jo ~ 120 for the current at the recon- 
nection layer central point O [isl [l6| . Our approxi- 
mate Eqs. ([2T|) . derived for the localized resistivity case, 
give jo ~ {Brn/lri){VAU-i/'noY^'^ = 146, which is close to 
jo ~ 120 observed in the simulation. The small disagree- 
ment could be due to plasma compressibility, or due to 
the finite Lundquist number used in the simulation and 
due to the approximate nature of our theoretical model. 

We conclude that our theoretical model for mag- 
netic reconnection is in an agreement with the simula- 
tions by Baty et al. [3]. Our model is in reasonable 
agreement with several other previous numerical simula- 
tions of reconnection with spatially nonuniform resistiv- 

ity 0,i[iil,[i3. 



IV. DISCUSSION 

Let us now discuss our major results. 

First, equation pT|) that determines the reconnection 
current jo and the quasi-stationary reconnection rate 770J0 
is derived by using the "local" equations theoretical ap- 
proach (see the last column of Table [5]) . This approach 
involves only "local" quantities and equations that are 
defined on the interval OM across the reconnection layer, 
refer to Fig. [1] Thus, the quasi-stationary reconnection 
rate in a thin two-dimensional layer is determined lo- 
cally, in the layer central point O and in the upstream 
region outside the layer at point M. In other words, 
the rate is fully determined by a particular functional 
form of anomalous resistivity ri(jz,x,y) and by the lo- 
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cal configuration of the magnetic field in the upstream 
region. The later determines field and its scale 
L= [-2B,„/(a^By)™]i/2 at point M. As a result, the 
global properties of the reconnection layer (e.g. its length 
L' , the plasma outflow velocity V^^^ and the presence or 
absence of Petschek shocks) do not directly matter for the 
quasi-stationary reconnection rate. This is a very impor- 
tant result because the reconnection layer global geome- 
try can be very complicated (e.g. in turbulent plasmas). 
Of course, there exists an indirect dependence of the re- 
connection rate on the field global properties because the 
local field configuration in the upstream region is deter- 
mined by the field global configuration. The above state- 
ments are also true for all other "local" parameters - for 
the layer thickness So and for the reconnection velocity 
Vr, refer to Eqs. ([I2D and ([T3]). 

Second, the reconnection rate can also be estimated 
by using the "global" equations theoretical approach, in 
which the whole reconnection layer is considered, includ- 
ing the downstream region at point O' (refer to Fig. [1]). 
The equations used in this approach are presented in the 
second column of Table [HI We would like to point out 
that two of these equation play the key role in correct de- 
termination of the reconnection layer length L' and geom- 
etry. The first key equation is the jump condition on the 
Petschek shocks, B'^/y/p « « Vr . whose importance 
was first pointed out by Kulsrud The second key 

equation is the constancy of the electric field z-component 
along the reconnection layer, rjojo = v'j'z ~^^out^xj which 
has been overlooked in the previous theoretical models 
(e.g. in the Petschek model). 

Third, for the case of a strong dependence of re- 
sistivity on the current, i.e. when {jo/'>h){dri/djz)o ^ 
max[l, L^/Z^], our results, given by Eqs. ([20|. coincide 
with the results obtained by Kulsrud [iJI- Thus, both 
the "local" and "global" theoretical approaches confirm 
Kulsrud's results and ideas, contrary to the doubts raised 
in paper by Baty et al. [15,]. 

Fourth, we found that in the case of uniform or quasi- 
uniform resistivity the magnetic reconnection rate is the 
slow Sweet-Parker rate and not the fast Petschek rate, 
see Eqs. ([T5|) . This theoretical result follows from rig- 
orous analytical derivations and agrees with numerical 
simulations. At the same time it contradicts the origi- 
nal Petschek theoretical model. Let us consider both the 
"global" and "local" analytical approaches in the case of 
a strictly uniform resistivity i] = const = 770, and let us 
explain why the Petschek reconnection layer geometry is 
not realized in this case. We take the "global" equations 
approach first. In the case of constant resistivity Eqs. ^ 
and pO)) result in equation (770^0 — r]'j'^)/r]ojo ~ {L'/L)^ 
for the fractional drop of the rjjz term along the recon- 
nection layer. On the other hand, the constancy of the 
electric field z-component along the reconnection layer 
implies that this fractional drop is {rjojo — v'j'z)/Vojo ~ 
V^^^B'^/rjojo ~ 1, where we use the "global" equations 
on lines 3 to 6 in the second column of Table [6l These 



two equations agree only if L' k, L, which means that 
the geometry of the reconnection layer is Sweet-Parker 
and not Petschek. Next, we take the "local" equations 
approach. In this "local" approach we prefer not to con- 
sider the reconnection layer geometry and any "global" 
parameters, such as the layer length L', for the calcu- 
lation of the reconnection rate. Instead we argue as 
follows. Refer to the "local" equations in the last col- 
umn of Table [6l Equations on lines 1, 2 and 5 result 
in an estimate of the plasma outflow velocity derivative 
{dyVy)o « Vojo/^nn see Eq. Plasma acceleration 

equation on line 3 results in the upper estimate for the 
Bx field derivative, {dyBx)o < p{dyVy)o^ / jo- As a result, 
we can find the upper estimate for the reconnection cur- 
rent jo from the "local" equation on line 6, which is the 
condition of constancy of the electric field z-component 
along the reconnection layer. If resistivity is uniform, 
this estimate turns out to be the Sweet-Parker value, 
Jo ~ {B,„/L){VaL/tJo)-^/\ as given by Eq. CID. In 
addition, if one estimates the reconnection layer length 
as L' ~ B'^/{dyBx)o ~ VA/{^yVy)o^ he/she would again 
recover the Sweet-Parker result L' ^ L |lj| . We conclude 
that in the case of constant resistivity both the "global" 
and "local" equations approaches consistently lead to the 
Sweet-Parker reconnection rate and the Sweet-Parker ge- 
ometry of the reconnection layer (L' « L). 

Finally, let us point out that whether reconnection is 
unforced (free) or forced does not matter for our results 
and conclusions. Indeed, on one hand, in the case of 
unforced reconnection one solves Eqs. pTjl and (fT5)) for 
the current jo and for the reconnection velocity Vr. The 
solution will depend on the magnetic field Bm in the up- 
stream region at point M. On the other hand, in the case 
of forced magnetic reconnection the reconnection velocity 
Vr is prescribed and fixed. In this case the field Bm in 
the upstream region should be treated as an unknown pa- 
rameter, and Eqs. (|lip and ([T^ are to be solved together 
in order to find the correct quasi-stationary values of jo 
and B„i. In other words, in the forced reconnection case 
an initially weak outside magnetic field Bm gets piled 
up to higher values until the resulting current jo in the 
reconnection layer becomes large enough to be able to 
match the prescribed velocity Vr of magnetic flux and 
energy supply in the upstream region. 
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